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DESIGN AND OVERALL PERFORMANCE OF A TWO-STAGE AXIAL-FLOW 
PUMP WITH A TANDEM-ROW INLET STAGE 
by Donald M. Sandercock and James E. Crouse 
Lewis Research Center 

SUMMARY 

The inlet stage group for a high suction specific speed, multistage axial-flow pump 
was designed and tested in water. Inlet flow vector diagrams were optimized for suction 
specific speed at a suction specific speed of 30 000. The first stage had a tandem ar- 
rangement of a cambered, tapered inducer and a transition rotor that adds enough head 
to keep the following highly loaded stages free of cavitation at design flow. The first- 
stage stator set up flow conditions suitable for the highly loaded rotor. 

Under noncavitating conditions, the pump performance at design flow was below the 
anticipated head and efficiency because the transition and second-stage rotors failed to 
achieve anticipated energy addition and efficiency. As the flow coefficient was reduced 
below 0. 9, the pump developed regions of reverse flow. The characteristic had little 
drop in head rise in the stall range. The off-design blade row characteristics showed 
major redistributions of work between rotors. 

At design flow and a suction specific speed of 27 000 in water, the pump head rise 
had begun to fall from the noncavitating values. 

INTRODUCTION 

The primary function of inlet stages of a high suction specific speed multistage axial- 
flow pump is to accept cavitation but continue to produce sufficient pressure to suppress 
cavity formation and prevent its occurrence in main pressure producing portions of the 
pump. The inlet stage experiences significant amounts of two-phase flow; its first blade 
row, or inducer, characteristically includes high blade stagger angles (as measured 
from the axial direction), a small number of blades, large solidities, low hub-tip radius 
ratios, and low blade loading. In contrast, the design of the high pressure producing 
stages with single-phase flow employ lower blade stagger angles (50° to 70°), lower 






solidities, larger numbers of blades, higher hub-tip radius ratios, and limiting values, of 
blade loading. These two types of stages are linked by means of one or more transition 
stages. Minimum pump weight will be realized when the number of transition stages is 
minimized. 

In this study, the usual inducer and a transition stage have been combined into a 
single stage composed of a tandem arrangement of two rotating blade rows followed by a 
single stator row. Thus the inducer rotor, the transition rotor, and the stator row com- 
prise the first stage. The test pump configuration is completed by the addition of a rotor 
of the type used in the high pressure producing portion of a multistage axial-flow pump. 
The first rotating blade row is typical of inducers including a low level of loading, a low 
hub-tip radius ratio, and several long- chord blades. The tandem- mounted transition 
rotor utilizes conventional blade shapes and chord lengths together with an increased level 
of loading. Such an arrangement would eliminate one stator row with the sacrifice of the 
cavitation reduction effects of the stator static pressure rise; alternately, as compared 
with the rear portion of a single extended chord rotor, the transition rotor starts its load- 
ing with clean surface boundary layers. 

This report presents design procedures and overall performance of this multistage 
pump under noncavitating and cavitating flow conditions. The variation of head rise 
across the midspan of each blade row as flow is reduced through the stall region to shut- 
off is also presented. 


PUMP DESIGN 

Overall Design Considerations 

At the initiation of this design, it was anticipated that the pump evolved would be 
operated in both water and liquid hydrogen. Consequently, operation in both fluids and 
the capabilities of water and liquid hydrogen test facilities were considered in selecting 
design parameters. 

It was desired to obtain a suction specific speed of 30 000 with this pump. From 
optimization procedures similar to those outlined in reference 1, a suction specific 
speed of 30 000, a hub-tip radius ratio of 0. 4, and the assumption of no inlet whirl com- 
bines to give an optimum relation between the axial velocity and blade speed of 


v z 

— = 0. 108 
U 


(1) 


(Symbols are defined in appendix A. ) If equation (1) is applied at the blade tip location, 
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tins relation becomes, by definition, the ideal inlet flow coefficient, or 
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The selection of an outer diameter of 6. 50 inches was based primarily on considerations 
of liquid hydrogen test facility capability and acceptable length of operating time in liquid 
hydrogen. Flow rates were calculated by combining the design flow coefficient and the 
hub-tip radius ratio of 0. 4 with a number of typical operating speeds in liquid hydrogen. 

It is desirable from both cavitation and blade loading considerations that the flow co- 
efficient increase (area decrease) when progressing from inlet stages with high suction 
specific speeds and low blade loading to latter stages with lower suction specific speeds 
and higher blade loading. In this design the average flow coefficient is increased by a 
factor of two from the inlet of the first rotating blade row (inducer) to the inlet of the 
second-stage rotor. Half of this increase is accomplished across the inducer and the re- 
maining across the transition rotor of the first stage. No change in average flow coef- 
ficient was used across the second-stage rotor. 

A preliminary sketch of the flow path contour based on the aforementioned area var- 
iation, a constant tip diameter, typical axial depths for each blade row, and sufficient 
spacing between blade rows to permit survey instrumentation to be used was constructed. 
The sketch showed a severe inducer hub curvature and a hub-tip radius ratio of 0. 761 
with a blade height of 0. 776 inch for the second-stage rotor. In order to alleviate the 
hub curvature (and its affects on radial equilibrium requirements) and to improve the 
blade height of the second and any succeeding stages, the outer wall casing was tapered 
linearly across the inducer row to 0. 95 of the inlet diameter (from 6. 500 to 6. 175 in.). 
For the same flow coefficient variation across the first stage, the hub-tip radius ratio of 
the second-stage rotor becomes 0. 720 with a blade height of 0. 865 inches. A sketch of 
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this flow passage showing axial location of measuring stations is presented in figure 1? 

In this design, the head rise coefficient of each rotating blade row was constant 
from hub to tip. Radially constant total pressure tends to minimize mixing losses and 
reduce radial gradients of axial velocity thus easing matching problems with succeeding 
blade rows. Additional guide lines used to select approximate design levels of blade 
loading or head rise coefficient, radial distributions of loss coefficient, and type of ve- 
locity diagram in the various blade rows were as follows: 

(1) A design head rise coefficient of approximately 0. 120 was selected for the in- 
ducer. This value is consistent with the level of performance generally obtained from 
inducers with good cavitation performance. This head rise coefficient would provide a 
cavitation number (see eq. (B12a)) of approximately 0. 3 at the inlet to the transition 
rotor tip element. 

(2) Blade loading across the transition rotor was limited to a blade diffusion factor 

of approximately 0. 4 in the blade tip region and 0. 6 at all other radial locations. Velocity 
diagrams at the transition rotor outlet were calculated from the inducer inlet velocity dia- 
grams and a design head rise coefficient and a radial distribution of loss coefficient as- 
sumed to occur across the two rotating blade rows. 

(3) Velocity diagrams for the second- stage rotor utilized some inlet whirl at all 
radial locations. This reduces the cavitation number requirements at the rotor inlet and 
tends to equalize the portion of overall stage diffusion accomplished across the rotor and 
stator blade rows. A rather high design blade diffusion factor of approximately 0. 55 at 
the blade tip was selected. 

Having determined pump geometry and a preliminary value of overall head rise co- 
efficient, a pump rotational speed of 5415 rpm for operation in water was selected. 
Principal factors affecting this selection were the limiting pressures of the collector 
downstream of the pump and the need for obtaining as high an inlet pressure correspond- 
ing to the 30 000 suction specific speed operating condition as possible. 


Detailed Blade-Row Design 

The pump detailed design approach is based on the assumptions of blade element 
flow and axial symmetry. The blade element approach assumes that flow in the circum- 
ferential plane may be described by considering the flow about blade profiles formed by 
the intersection of a flow surface of revolution with pump blading. Axial symmetry as- 
sumes that the flow conditions in the circumferential, or blade-to-blade, plane can be 
represented by an average value. Equations describing the radial variations of these 
values are written from continuity, energy addition, and pressure equilibrium from hub 
to tip. 
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* In this design, the blade elements were assumed to lie on conical surfaces intersect- 
ing the blade inlet and outlet computational planes at radii enclosing equal portions of 
weight flow. Simple radial equilibrium (refs. 2 and 3), as defined by 

y2 

^ ( 2 ) 

3r gr 

was applied in all design computations. The form equation (2) takes when used in a de- 
sign application is given in appendix B together with other design equations. 

Velocity diagrams were computed at the inlet and outlet of each blade row at 11 rad- 
ial locations that divided the total flow into 10 equal portions. The large number of flow 
surfaces used reflects the convenience of being able to program the design equations for 
solution on a high-speed computing machine rather than an accuracy requirement. Design 
parameter together with velocity diagram calculations based on a pump inlet tip speed of 
153. 5 feet per second are tabulated in tables I (p. 6) and II (p. 7) for a selected number 
of blade elements. 

At the first rotor inlet the velocity diagrams are based on the assumptions of (1) no 
inlet whirl (V^ j = 0) and (2) a radially constant inlet flow coefficient. 

The radial distribution of inducer outlet velocity diagrams was calculated from the 
equations listed in appendix B. Known inputs to the equations, include flow rate, passage 
hub and tip radii, inlet velocity diagrams, and selected radial distributions of head rise 
and loss coefficient. The latter were selected primarily on the basis of loss distributions 
measured during operation of a low hub-tip radius ratio rotor in water. 

From the inducer blade trailing edge (inducer outlet computational plane) to the 
transition rotor blade leading edge (transition rotor inlet computational plane), a signifi- 
cant change in flow area occurs due to hub slope and the relatively large axial clearance 
space. In order to obtain more accurate inlet flow conditions to the transition rotor, 
the inducer outlet velocity diagrams were corrected by accounting for change in area 
and conservation of angular momentum with change of radius. The velocity diagrams 
at the transition rotor outlet were computed from the inducer inlet diagrams, the desired 
head rise to be produced by the inducer and transition rotors, and an assumed radial dis- 
tribution of loss coefficient that would occur in the flow across the two rotating rows. 

This overall loss coefficient distribution was calculated from the loss coefficient distri- 
butions assumed for the inducer and transition rotors individually. The latter were 
chosen from experience with rotors operating in water (e. g. , ref. 4). 

The first- stage stator must accept the outlet flow from the transition rotor and turn 
it to the direction required by inlet velocity diagrams of the second-stage rotor design. 
Hence, the stator design is dependent upon the designs of both the preceding and succeed- 
ing rotors. 

The conditions chosen for the second-stage rotor design were as follows: 
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TABLE n. - BLADE ELEMENT DESIGN PARAMETERS 


Blade 

row 

Stream- 

line 

number 

Radius, 

r, in. 

Entrance 

flow 

coefficient, 

[U t =153.5] 

Rotor head 
rise coeffi- 
cient, 

* 

(u t = 145.8] 

Rotor 

effi- 

ciency, 

V 

Blade 

diffu- 

sion 

factor, 

D 

Loss 

coeffi- 

cient, 

“CO 

Relative 
change 
in flow 
angle, 
Aj3* 

Inlet 

Outlet 

1R 

0 

3.250 

3.087 

0. 108 

0. 133 

0. 655 

0.2635 

0. 125 

3.2 


2 

2. 964 

2. 844 





.740 

.2622 

. 100 

6.0 


4 

2. 647 

2.626 





. 809 

.2496 

.084 

8. 0 


6 

2.287 

2.411 





.856 

.2309 

.080 

9.8 


8 

1. 859 

2. 189 





.899 

. 1871 

.080 

11.6 


10 

1.296 

1. 954 

1 


L_ _ 


. 940 

. 0568 

. 090 

12. 6 

2R 

0 

3. 087 

3.087 

0. 155 

0. 227 

0. 896 

0.4171 

0.079 

2.3 


2 

2. 890 

2. 864 

. 181 



. 933 

.4356 

.055 

9.4 


4 

2.705 

2.700 

. 196 



. 948 

.4423 

.049 

14.6 


6 

2. 520 

2. 550 

.205 



. 956 

.4511 

.048 

19.0 


8 

2.329 

2.404 

.212 



. 960 

.4564 

.052 

23. 8 


10 

2. 130 

2.258 

.220 



.963 

.4498 

.060 

29.4 

IS 

0 

3.087 

3.087 

0. 128 

— 

__ 


0. 5884 

0. 107 

41. 8 


2 

2. 864 

2.950 

. 198 

--- 

-- 


. 5757 

.096 

28. 5 


4 

2.700 

2.800 

.237 

--- 

-- 


. 5681 

.088 

20.7 


6 

2. 550 

2. 635 

.259 

--- 

-- 


. 5613 

.080 

16.0 


8 

2.404 

2.456 

.277 

--- 

-- 


. 5520 

.074 

12.3 


10 

2.258 

2.258 

.293 

— 

— 


. 5403 

.067 

9.7 

3R 

0 

3.087 

3.087 

0.216 

0.340 

0. 887 

0. 5573 

0.113 

9.9 


2 

2. 864 

2. 940 





.914 

. 5692 

.082 

12.9 


4 

2.700 

2.785 





. 933 

.6160 

.063 

17. 9 


6 

2. 550 

2.622 





.949 

.6822 

.053 

26.3 


8 

2.404 

2.447 





. 962 

.7532 

.051 

41.3 


10 

2.258 

2.258 





. 972 

.7853 

.070 

66. 0 


(1) Radially constant head rise 

(2) Radially constant axial velocity 

Satisfying the aforementioned conditions means that a symmetrical velocity diagram can 
be attained across only one blade element. In this design, the symmetrical velocity 
diagram was located at the mean blade element. The blade diffusion factor was limited 
to a value of approximately 0. 55 at the tip but no limiting values were applied to other 
blade elements. The blade loading over the complete blade span was higher than that 
used in current practice. The selected radial distribution of loss coefficient was based 
on preliminary performance data from similar blade rows but with lighter loading. Using 
the previous conditions together with the definitions of a symmetrical velocity diagram 
results in the following velocity diagram calculations: 
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(1) At the blade mean streamline, a symmetrical velocity diagram is computed (see 
ref. 2 for discussion of symmetrical velocity diagram) based on an assumed blade diffu- 
sion factor. From the assumed loss coefficient, a head rise is calculated for application 
at all elements. 

(2) The radial equilibrium equation (eq. (B5c)) is used to calculate inlet values of 
tangential velocity at all other radii. 

Substitution of aforementioned design conditions (9 AH/9r = 0, 9V / 9r = 0, and Zo = known 

z 

value) into equation (B5c) leaves rotor inlet (stator outlet) tangential velocity as the only 
unknown if blade element flow is assumed. Several iterations of the aforementioned 
procedure may be necessary before an acceptable value of blade tip diffusion factor is 
obtained. Iterations are also generally necessary in steps (1) and (2) before the calcu- 
lated radial distribution of axial velocity at the stator outlet approaches the assumptions 
of axial velocity necessary in the computational procedures. 

Throughout the design procedures used herein, the radial distributions of loss co- 
efficient have been assumed. In more recent design procedures, the loss coefficient 
enters the calculations as some function of the blade loading and radial location of the 
blade element. 

After the velocity diagram design data had been completed and applied in blade de- 
sign calculations, a discrepancy in the second-stage rotor inlet velocity diagram was 
found. Across the stator, the assumed radial distribution of loss coefficient and the 
radial gradients of tangential velocity as established by the transition rotor outlet and 
second-stage rotor inlet flow conditions result in some mass flow shifts. Thus blade 
elements across the stator intercept inlet and outlet computational planes at different 
radii. The velocity diagram information computed for stator outlet (station 4 radii) was 
applied to the succeeding rotor inlet velocity diagram calculations at the stator inlet 
radii (station 3) rather than at the correct stator outlet values (station 4). This is re- 
flected in the blade speeds presented for axial station 4 in table I (p. 6), which correlate 
with the radii at axial station 3 rather than 4. Subsequent tables of design values also 
show differences in the location of the streamlines leaving the stator outlet and entering 
the succeeding rotor. This results in maximum errors of 0. 9° in relative flow angle and 
5 feet per second in relative velocity, which occur at approximately the blade midspan. 
Because fabrication of the pump was nearly complete and the degree of error was small, 
it was not deemed necessary to correct these small discontinuities in streamline paths. 


Blade Design 

Blade shape. - Originally the inducer utilized a circular arc meanline with constant 
thickness from inlet to outlet along a given radius. The blade thickness varied linearly 
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with radius, increasing toward the hub. The circular arc meanlines were assumed to 
lie on conical surfaces intersecting the blade inlet and outlet computational planes at 
radii dividing these planes into equal flow areas, that is, blade element surfaces. 

Before the inducer was built, some unpublished data were obtained on a similar in- 
ducer indicating that the performance of this inducer could be improved by revising the 
axial distribution of flow area to increase the passage area in the inlet portions of the 
blade row. This modification was expedited by adding a helical blade section to the 
leading edge of the inducer rotor. The helical portion extended an axial distance of 
0. 8 inch upstream of the inducer leading edge and had a constant passage area. At the 
blade tip (r = 3. 25 in. ) the helix angle matched the rotor blade inlet angle. At other radii 
some degree of fairing from the helical portion to the initial rotor design was required. 

All other blade rows employed double circular arc blades. In all radial locations 
the double circular arcs lay along the assumed blade element, or streamline, surfaces. 

Significant blade design parameters for six blade elements are presented in table III 

(p. 10). 

Incidence and deviation angles. - For this design, selected incidence angles were 
based on some limited data from helical inducers and a double circular arc rotor. 

Deviation angles of all the double circular arc blades rows were computed from the 
deviation angle rule of reference 5, which is defined as 


5 = 



( 3 ) 


where m for circular arc cambered airfoils is given as a function of blade setting angle 
up to a maximum value of 60°. Because the blade shape and range of inlet blade angles 
of the inducer were very different from those from which the design rule was formulated, 
the rule was not applied; deviation angles equal to one-fourth of the required turning 
angle were assigned. The latter values compare favorably with values computed from 
equation (3) if m is extrapolated to the inducer blade setting angles. It also can be 
noted from table 131 that the setting angles of a number of blade elements with double cir- 
cular arc profiles were somewhat greater than the maximum values covered in the for- 
mulation of the rule and extrapolation was necessary. 

Blade coordinates. - In the design procedures discussed thus far, each blade element 
has been assumed to lie on a conical surface intersecting the blade row inlet and outlet 
computational stations at points of equal weight flows. The complete blade is formed by 
stacking the blade elements on a radial line passing through the center of gravity of the 
individual sections. Cylindrical coordinates of the blade pressure and suction surfaces 
are obtained from a high-speed computer program. 

For fabrication purposes, it is desirable to provide blade surface coordinates along 
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TABLE HI. - BLADE DESIGN PARAMETERS 


row line 

number 


Radius, 

r, in. 

Inlet 

Outlet „ 

c 

3.250 

3. 087 

2.964 

2. 844 

2.647 

2.626 

2.287 

2.411 

1.859 

2. 189 

1.296 

1.954 

3.087 

3. 087 

2. 890 

2. 864 

2.705 

2.700 

2.520 

2. 550 

2.329 

2.404 

2.130 

2.258 

3.087 

3.087 

2.864 

2.950 

2.700 

2. 800 

2. 550 

2. 635 

2.404 

2.456 

2.258 

2.258 

3. 087 

3. 087 

2.864 

2. 940 

2.700 

2.785 

2. 550 

2. 622 

2.404 

2.447 

2. 258 

2.258 


Inci- 

Devia- 

Blade 

Blade 

Blade 

Blade 

Maxi- 

En- < 

Outlet 

Num- 

dence 

tion 

so- 

set- 

cam- 

chord, 

mum 

trance 

blade 

ber 

angle, 

angle, 

lidity, 

ting 

ber 

c, 

blade 

blade 

angle, 

of 

i> 

5, 

a 

angle, 

angle, 

in. 

thick- 

angle, 

K o> 

blades 

deg 

deg 


y, 

(p°> 


ness, 

K e’ 

deg 





deg 

deg 


^max’ 

deg 









in. 










(a) 




3. 0 

0. 8 

3.65 

— 

1.0 

— 

0. 060 

80. 9 

79. 9 

3 


0 4.7 

5.4 
6.7 
b 8.7 

.7 11.0 

4.0 13.0 


3. 17 

2.95 

2.82 

2.71 

2.62 

1.47 76.1 
1.58 68.2 
1.68 62.3 
1.79 56.9 
1. 92 51.0 
2.07 43.9 

1. 62 45. 9 
1.73 46.7 
1. 82 48. 0 
1. 93 49. 7 
2.06 51.5 
2.22 52.9 

1. 53 66. 5 
1.63 63.5 
1.73 59.2 
1.83 52.5 

1.95 41.9 
2.10 24.0 


? 74. 4 

.090 68. 1 

1.50 0.075 


55.4 
37.6 
27.2 
20.8 
15.9 
12. 5 


♦ . 120 
1. 50 0. 120 


♦ .060 
1.75 0.070 




a Linear distributions 
b r„ = m = 0.010 in. 


in blade maximum thickness used. 


s 



i 



Figure 2. - Two-stage axial-flow pump. 


/-Circulating pump for 



a constant radius plane. These are obtained 
using the blade surface coordinates from the 
computer program in a graphic procedure. 

A photograph of the pump is shown in fig- 
ure 2. 

APPARATUS AND PROCEDURE 
Test Facility 

The pump was tested in the Lewis water 
tunnel. The principal features of the water 
tunnel are shown in figure 3 and discussed 
in reference 6. Prior to the water tests, 
solid material larger than 5 microns was re- 
moved by the filtering system and the gas 
content was reduced to less than 1 part per 
million by weight with the auxiliary degas i- 
fying loop. During the tests, the gas con- 
tent was maintained below 3 parts per mil- 
lion. 
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55432 

(a) Total pressure claw. (b) static pressure wedge. (c) Total pressure wake rake. 

Figure 4. - Survey probes. 

Instrumentation 

The instrumentation used to obtain overall test loop conditions consisted of a venturi 
flowmeter, an electronic speed counter used in conjunction with a magnetic pickup, an 
automatic water temperature control and recorder, and a pressure transducer with the 
associated power and recording equipment to measure the loop pressure. The survey 
instruments (fig. 4) consisted of claw-type probes for measuring total head and angle be- 
hind rotating blade rows, a total head rake for measuring the wake behind the stators 
and wedge-type probes for measuring static head and angle. The measuring station lo- 
cations are indicated in figure 1 (p. 3), which shows a meridional section of the pump. 
Each claw and wedge had null -balancing stream- direction- sensitive equipment that auto- 
matically alined the probe to the direction of flow. A head calibration factor for each 
wedge static probe was determined in a low- speed air tunnel at approximately the same 
Reynolds number used in water and applied to the static heads measured in the water 
tunnel. The inherent accuracies of measurement and recording devices were as follows: 


Flow rate, Q, percent <±1. 0 

Rotative speed, N, percent ±0. 5 

Differential heads, AH, percent ±1.0 

Flow angle, |3, deg ±0. 5 


This does not include possible errors due to circumferential variations in flow, unsteady 
flows, and other time or space gradients. 
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- TABLE IV. - OPERATING 
CONDITIONS 


Procedure 


Rotor tip speed, 
u f 

ft/sec 

Net positive 
suction head, 

H sv- 

ft 

123 

246 

154 

248 

154 

135 

154 

38 

154 

15 


Characteristic performance curves were obtained by 
operating the rotor at constant rotative speed and net posi- 
tive suction head over a range of flow control valve settings. 
The data were taken by surveying total and static heads, and 
flow angles across the annulus at preprogramed radial posi- 
tions located at approximately 10, 30, 50, 70, and 90 per- 
cent of design flow from the tip. The data presented were 
taken along characteristic curves for the conditions listed 
in table IV. The equations used to calculate selected per- 


formance parameters are presented in appendix B. 


Data Reliability 

Some measure of data reliability is available from weight flow checks, which are 
defined as the difference between the flow calculated from the survey data and the venturi 
flowmeter reading divided by the venturi reading. The weight flow checks are plotted 
against pump inlet flow coefficient in figure 5 for each of the five measuring stations. 

The circles in figure 5 (p. 14) represent data taken under noncavitating conditions and 
other symbols represent data taken under various degrees of cavitation. 

The noncavitating data weight flow checks are within ±6 percent at measuring stations 
1, 2, 3, and 5 with the exception of a few low-flow points where there is a flow reversal 
at the inlet of the inducer tip. At the measuring stations behind rotors, good blade-to- 
blade averages are obtained because the blades are moving with respect to the instru- 
ment. At station 4 behind the stators the weight flow differences run as high as 14 per- 
cent. At these measuring stations, the wakes are not moving with respect to the instru- 
ments so good average conditions are somewhat difficult to obtain. With the total pres- 
sure rake a good average value of total head is expected, but the blockage effect of the 
rake and particularly the angle measurements with the static probe present some prob- 
lems. Although better data accuracy is desired, the accuracy obtained is as good as 
can be expected in small multistage machines with normal test procedure. 

In general, the cavitating data weight flow checks are poorer than those calculated 
from the noncavitating data. One cause of this is probe cavitation, particularly on the 
static probe under some conditions. Observation of such cavitation through transparent 
casings over other rotors indicate that the total pressure readings are likely to be un- 
affected until the inlet pressure is lowered somewhat from that at which cavitation is 
first observed on the static probes. Thus, it is expected that the head rise measure- 
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ments under cavitating conditions are some- 
what better than the weight flow data indi- 
cates. The accuracy of the velocity and ef- 
ficiency calculations, however, can be ex- 
pected to deteriorate with cavitation. 


RESULTS AND DISCUSSION 

The test results are presented as non- 
cavitating and cavitating data. In this in- 
vestigation, the noncavitating flow regime at 
a given wheel speed is defined as operation 
at inlet pressures above the pressure where 
a change in overall performance of the pump 
with a change in inlet pressure cannot be 
detected. Thus, the pump is not considered 
to have entered into the cavitating regime 
until cavitation is extensive enough to affect 
overall performance. 


Noncavitating Performance 

Overall pump performance. - The over- 


Figure 6. - Overall pump performance for noncavitating 
conditions. Net positive suction head, 246 feet. 


all noncavitating performance is shown in 
figure 6 as head rise coefficient and effi- 
ciency plotted as functions of inlet flow co- 
efficient. Design flow coefficient is considered to be 0. 104. This value is lowered from 
the ideal flow coefficient of 0. 108 to allow for hub and tip boundary layer blockage. The 
head rise coefficient and efficiency represent mass-averaged values, which were com- 
puted by the methods of appendix B, whereas the flow coefficient is an average value 
based on venturi- measured flows and inlet geometric area. The data shown were taken 
at two rotative speeds, 4330 and 5415 rpm, representing blade tip speeds of 123 and 
154 feet per second, respectively. The data at 4330 rpm represent the noncavitating 
condition as defined earlier over the whole range of operation. Only the lower pressure 
(high flow) end of the characteristic curve at the design speed of 5415 rpm is shown. At 
low flow coefficients, the high inlet pressure combined with the head rise of the pump 
raised the exit pressure level to the point where the collector spread enough to open a 
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face seal from its mating surface. Even at this relatively high inlet pressure a falloff k 
in performance due to cavitation is indicated at high flows, but the 5415-rpm data fair in 
with 4330-rpm data near the maximum flow point and are presumed to follow the non- 
cavitating characteristic over the rest of the range. The 5415-rpm data are shown on the 
noncavitating curve because it is the design speed and it is used along with the 4330-rpm 
data as a reference condition for the lower inlet pressure data at 5415 rpm. At reduced 
inlet pressures, the collector pressure level was low enough so that the pump could be 
safely operated at the higher speed over the complete flow range. 

The range of operation shown in figure 5 is limited on the high-flow end by the test 
loop losses when the flow control valve is completely open. As the flow was reduced, a 
recirculation at the inlet tip of the inducer began at a flow coefficient of 0. 09 and in- 
creased in radial depth as the flow was decreased to 0. 07 where radial surveys were 
terminated. The recirculation is probably due to a radial equilibrium requirement rather 
than blade stall because it was not accompanied by any appreciable increase in rig vibra- 
tion. Since the presence of a recirculating flow region of some magnitude has the poten- 
tial of triggering instabilities in some systems, however, the point of initiation of re- 
circulation (cp = 0. 09) was selected as the criterion for the limit of the stable operating 
range. When recirculation at a measuring station occurs, mass-averaged values of head 
coefficients and efficiencies are difficult to determine. Head coefficients and efficiencies 
are mass weighted on the blade row outlet axial velocity distribution, but difficulties are 
still encountered when a recirculation occurs at the inlet only, because the inlet probes 
measure work done on the fluid by the blade row in the recirculation region. In general, 
the head values obtained from mass averaging under such adverse conditions are quite 
close to the mid survey point head rise values but the mass-averaged values of efficien- 
cies are doubtful. It is in this light that the mass-averaged values for a flow coefficient 
below 0. 09 are presented. More information on the characteristic curve in the low-flow 
region will be presented in a later figure in which midpoint data to shutoff flow is shown. 

In comparing the data of figure 6 with design, it is apparent that the design head rise 
and efficiency were not attained. Not only was the efficiency low but the energy addition 
was not achieved since the percentage deficiency in head rise is somewhat larger than 
that indicated by the decrement in efficiency. Discussion of this is included in the blade- 
row performance section. Another point of interest is the location of design flow on the 
characteristic curve with respect to efficiency and range. Since the efficiency curve is 
relatively flat in the region of design flow, the design point could be moved to a slightly 
higher flow coefficient without a sacrifice in efficiency if it were desirable to have addi- 
tional flow margin from design to the low flow limit of the stable operating range. 

Blade-row performance . - The blade-row performance is presented in figure 7 by 
plotting cumulative head rise through the blade rows and efficiencies of the individual ro- 
tating blade rows against pump inlet flow coefficient. At design flow, the inducer 
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(a) Cumulative head rise. (b) Rotor efficiency. 

Figure 7. - Noncavitating performance. Net positive suction head, 246 feet. 

achieved its design point and the stator losses were low; however, the transition and 
second-stage rotors fell far short of the design goals. The head rise across the transi- 
tion rotor was 0. 828 of the design value, and the efficiency was 0. 945 of the design value. 
From these figures, the energy addition is only 0. 876 of the design value, which is the 
major reason the transition rotor failed to achieve design head rise. The second- stage 
rotor head rise was 0. 760 of the design value, and the efficiency was 0. 860 of the design 
value. Thus, the energy addition was 0. 883 of the design value, and hence, the energy 
addition and efficiency defects are about equally responsible for the head deficiency of 
this rotor. 

From overall data, it is difficult to clearly resolve the causes of the deficiencies 
indicated above; however, some general observations are made. In regard to the transi- 
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tion rotor, it is suspected that the energy addition was not achieved because the deviation 
angles were larger than those anticipated in the design. It appears unlikely that flow 
separation on blade surfaces was the major cause since the blade loading was not espe- 
cially high and the efficiency did not decrease as significantly as the energy addition. A 
more likely reason is associated with the fact that, as was noted previously in the design 
section, blade setting angles of a number of blade elements (in tip region) were larger 
than the values used to formulate the deviation angle prediction rule used in the design. 

Across the second stage rotor with its high design blade loading, the significantly 
lower than design measured efficiency indicates the existence of thickened blade surface 
boundary layers or possibly some flow separation. Attendant effects would be to increase 
deviation angles. In addition, the blade setting angles for blade elements in the tip region 
required some extrapolation of the deviation angle design rule. Thus, although it is be- 
lieved that the second-stage rotor did not achieve the design energy addition because it 
did not produce the desired fluid turning, the extent to which the deficit should be attrib- 
uted to flow separation and extrapolation of the design deviation angle rule could not be 
further evaluated. 

The blade-row characteristics in figure 7 show how the distribution of work between 
rotors changes with flow coefficient. At high flows the transition rotor does practically 
all the work and the inducer is turbining. A zero inducer head rise was measured at flow 
coefficients greater than 0. 14 because the instrumentation was not set up to measure 
pressure drop. Since the transition rotor head rise and stator head drop are measured 
as differentials and added to the inducer head rise, the head levels in these regions are 
shown by the half -filled symbols in the figure. These were determined by assuming the 
inducer retains the same slope characteristic as it turbines. The stator losses in this 
region are several times the losses at design flow, and the second-stage rotor head rise 
barely compensates for these stator losses. As the flow is reduced, the inducer and 
second- stage rotor rapidly pick up large percentages of total pump work, while the tran- 
sition rotor head coefficient changes very little. With no stator between the inducer and 
the transition rotor, the inlet relative fluid angle, and hence, incidence angle of the 
transition rotor changes very little with flow. Thus, the fluid turning and energy addi- 
tion also change very little over the range of operation. 

At the low flows (c p < 0. 09) the inducer develops a recirculation at the inlet tip. As 
mentioned earlier, mass-averaged values are difficult to determine in a meaningful way 
under such conditions. In general, head distributions do not vary widely from hub to 
tip, so the method of averaging or flow reversals do not greatly affect the value attained. 
The value of efficiency, however, is quite sensitive to the axial velocity distribution and 
averaging method because the energy addition changes a great deal radially. No further 
effort was made to determine efficiencies of points where recirculation occurs, since, 
for the most part, the recirculation is inefficient and operation under such conditions is 
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to be avoided as much as possible. In this 
low-flow region the other rotor efficiencies 
are on the decline but not because of recircu- 
lations. The stator losses increase rapidly 
as the stator incidence angles are increased 
far above design values. 

Cavitating Performance 

Overall pump performance . - The over- 
all pump performance map for a range of net 
positive suction heads is shown in figure 8 as 
plots of head coefficient and efficiency against 
flow coefficient. Cavitation affects perform- 
ance on the high-flow end of the characteristic 
at relatively high net positive suction heads. 

In fact, there is a slight head decrease at an 
H gv of 248 feet when the rotor tip speed is 
154 feet per second. This corresponds to a 
suction specific speed of approximately 3800. 
As the net positive suction head is reduced, 
the head and efficiency dropoffs become 
inlet flow coefficient, j-j greater and the maximum flow becomes less. 

Figure 8. - Overall performance of axial-flow pump under In each case however until an H of 
cavitating conditions. ’ ’ SV 

15 feet is approached, the characteristic re- 
mains unaffected for design and lower flows. At an H gv of 15 feet (suction specific 
speed of approximately 27 000), the maximum flow rate is only slightly above the design 
flow, and the head, but apparently not the efficiency, has begun to fall at design and 
lower flows. The performance obtained indicates that a suitable operating suction spe- 
cific speed ought to be a little lower than 27 000. 

At maximum flow with an H gy of 15 feet, the characteristic curve for the pump is 
practically vertical (no measurable change in flow) for head rise coefficients ranging 
from 0. 17 to 0. 45. At the lowest head rise (open throttle operating point) visual obser- 
vations made during operation with a transparent casing showed that the stators and 
second-stage rotor were cavitating more severely than the inducer. These blade rows 
were very likely operating in the cavitation breakdown region where small changes in 
inlet pressure result in significant changes in head rise. As the throttle was adjusted to 
increase back pressure, the head rise across the inducer and transition rotor increased 
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Cumulative head rise 
coefficient from inlet to - 

o Inducer outlet 

- □ Transition rotor outlet 

o First-stage outlet 

v Pump outlet 

Plain symbols denote noncavitating 

- conditions 

Tailed symbols denote U* * 154 ft/sec 
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(a) Cumulative head rise. (b) Rotor efficiencies. 

Figure 9. - Comparison of noncavitating and cavitating performance. 


at only a moderate rate, but the effect was sufficient to give a rapid reduction in cavita- 
tion in the stator and second-stage rotor blade rows. The second-stage rotor head rise, 
and consequently the overall head rise, increased significantly without a measurable de- 
crease in inlet flow coefficient. The mechanism described is also validated by visual 
observations. During operation at an overall head rise coefficient of 0. 45, very little 
cavitation was observed in the stator or second-stage rotor blade passages. 

Blade-row performance. - The cumulative head coefficients and individual rotor ef- 
ficiencies for operation at an H of 15 feet are compared with the noncavitating data in 

bv 
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figure 9. At this H and design flow, a drop in head occurs only in the inducer as ex- 
pected. All other cumulative head rises are displaced downward by approximately the 
amount of the inducer head defect. At design flow the inducer produces enough head to 
make the rest of the pump run virtually cavitation free. The rotor efficiencies at design 
flow, however, do not show entirely expected trends. The inducer efficiency drops 
7 percentage points, and the transition rotor efficiency drops slightly. The second- 
stage rotor efficiency appears to rise about 5 percentage points and to compensate for the 
lower inducer efficiency to give essentially the same overall pump efficiency achieved 
under noncavitating conditions. This apparent rise in efficiency, however, is not be- 
lieved to occur since the weight flow checks, especially those behind the stator under 
cavitating conditions, indicate data inconsistancies. 



Figure 10. - Pump midpoint characteristic to zero flow under noncavitating conditions. 
Rotative speed, 4340 rpm. 
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Operation to Shutoff 


* 


Midpoint data were taken to determine the general pump characteristic to shutoff 
since this type of data can be useful for analysis of starting problems. For this data the 
total pressure probes were set at the 50-percent streamline for design conditions and 
left there for the run. The probes were alined to the flow with the automatic direction 
sensing equipment used in the surveys. The total pressure rake behind the stators was 
not used. The data making up this midpoint, noncavitating characteristic are presented 
in figure 10. Over the same flow range, midpoint data show very nearly the same char- 
acteristic as the mass-averaged data in figure 7 (p. 17). As flow is decreased into the 
stall region only a slight dip on head rise coefficient occurs before a rather steep nega- 
tive slope develops to shutoff. Further, as flow is increased from shutoff, this inlet 
stage group unstalls smoothly along the same characteristic with no hysteresis. 

Some complications arose in the determination of the work division between the in- 
ducer and transition rotors at low flows. Angle and head measurements between the in- 
ducer and transition rotor indicate that there is a backflow extending to the midstream- 
line position. Such a condition makes the inducer head rise appear to be very high and 
the transition rotor appear to turbine. Since the instrumentation was not set up to meas- 
ure pressure drops across the transition rotor, the transition rotor head rise and the 
combined inducer and transition rotor data is missing for low flows. No attempt was 
made to determine the nature or extent of this recirculation by surveying the passage. 
The recirculation at the inducer inlet observed in the radial surveys did not extend to the 
midstreamline. Thus, the pump overall midpoint head rise, which was measured di- 
rectly as a differential between outlet and inlet heads, was not affected by the afore- 
mentioned difficulties. 


SUMMARY OF RESULTS 

An inlet stage group for a high suction specific speed, multistage axial-flow pump 
was designed and tested in water. Analysis of the overall performance under both non- 
cavitating and cavitating conditions and limited performance measurements with the 
pump operating in unstable flow regimes indicated the following principal results: 

1. At design flow (average inlet flow coefficient Tp ^ = 0. 104) under noncavitating 
conditions, the pump performance was far below the design head and efficiency. The 
inducer and stator operated to design predictions but the transition rotor and second- 
stage rotor did not. 

2. The blade-row characteristics show major redistribution of work between rotors 
with change in flow coefficient. 
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* 3. At design flow under cavitating conditions (suction specific speed of 27 000) the 

pump head rise had begun to fall from the noncavitating values but the efficiency re- 
mained virtually unchanged. The drop in head occurred in the inducer. 

4. At low net positive suction heads the pump characteristic on the high-flow side 
becomes almost vertical. At the low- head portion of the characteristic, the stators and 
second-stage rotor cavitate more severely than the inducer. 

5. Noncavitating midpoint performance data in the very low-flow range show that the 
pump head rise has little dropoff in the stall flow range. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 11, 1965. 
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APPENDIX A 
SYMBOLS 


c blade chord, in. 

D blade diffusion factor, eq. (B3) 

g acceleration due to gravity, 

32. 17 ft/sec 2 

H total head, ft 

AH blade-element head rise, ft 

H_„ net positive suction head, 

S V 

H e - V ft 

h static head, ft 

hy vapor pressure, ft 

i incidence angle (angle between 

inlet flow direction and tangent 
to blade mean camber line at 
leading edge) deg 

k cavitation number, eq. (B12a) 

N rotative speed, rpm 

Q flow rate, gal/min 

r radius, ft 

t blade thickness, in. 

U rotor tangential velocity, ft/sec 

V absolute velocity, ft/sec 

/3 flow angle with respect to axial 

direction, deg 

A/3 change in flow angle 

y blade setting angle (angle between 

blade chord and axial direc- 
tion), deg 


6 deviation angle (angle between 

outlet flow direction and tangent 
to blade mean line at trailing 
edge), deg 

Tj rotor efficiency, percent 

k blade angle (angle between tangent 

to blade mean camber line and 
axial direction), deg 

0 blade solidity (ratio of blade chord 

to blade tangential spacing) 

<p flow coefficient, V z / (U t ) 

<p° blade camber angle, deg 

xp rotor head rise coefficient, 

g AH/(U 2 ) 

1 o 

co loss coefficient, eq. (B4a) 

Subscripts: 

e entrance to blade row 

h hub 

1 ideal 

He blade leading edge 

max maximum 

n axial station, see fig. 1 

o outlet of blade row 

t tip 

te blade trailing edge 

v measured with venturi meter 

z axial component 
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* 

6 tangential component 

1 inducer inlet measuring station 

2 inducer outlet and transition rotor 

inlet measuring station 

3 transition rotor outlet and stator 

inlet measuring station 


4 stator outlet and second- stage 

rotor inlet measuring station 

5 second-stage rotor outlet measuring 

station 

Superscripts: 

averaged quantity 
* relative 


V 
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APPENDIX B 


EQUATIONS 

Blade Element Equations 


Ideal head rise (energy input): 


U V 


AH i = 


o 0,0 




g 


( y ;> 2 

2g 


V’ 


(Bl) 


Efficiency: 


V = 


AH 


AH. 


(B2) 


Diffusion factor: 


p_, , r o V 0 , o - r e V e . e , t v ;, r 0 v e, 0 - r e v e,e 


v ; 


2cv;r 


V’ (r + r )ctV* 

e ' o e' e 


Total head loss coefficient: 




- H’ 


Kf 

2g 


When applied to a rotating blade row, equation (B4a) becomes 


AH. - AH 

U> = 

(Kf 

2g 

When applied to a stationary blade row, equation (B4a) becomes 


(B3) 


(B4a) 


(B4b) 
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CO = 


(B4c) 


H 



o 


In both design and analysis procedures, radial equilibrium was assumed to be 
achieved at each measuring station when the radial pressure gradient balanced the cen- 
trifugal acceleration of the fluid due to its whirling motion, or 



3r gr 


(B5a) 


Since 


h = 



2g 2g 




equation (B5a) can be written in terms of total head and velocities. With the radial 
component of velocity neglected, equation (B5a) becomes 




3V| 

Sr 



(B5b) 


For some applications, it is desirable to relate the radial distribution of outlet flow 
conditions to the inlet flow conditions and the change occurring across the blade row. By 
using the following relations: 


H = H + AH 
o e 


(vj 2 U v fl „ - u v a a (v ’) 2 

AH = AH i - co — — = t-l± * - co — — 

2g g 2g 


= i («o - U o V s, o tan Po ■ U e + U e v 2 , e tan »'e) ’ “ 


og; 

2g 


equation (B5a) can be written with a large variety of inputs. In the design procedure 
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used herein, equation (B5a) was written with AH, co , and inlet velocity diagrams as the 
known inputs. Hence, equation (B5a) takes the form 
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3H f 

o 0r o 
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(B5c) 

Later design procedures used AH^ as input and calculate a blade diffusion factor and an 
associated u>. 


Averaging Methods to Obtain Overall Parameters 

Allowance is made for the hub and casing boundary layers in the two end terms of 
each of the following summations: 


Mass-averaged head rise: 


r t t 

Z rV HAr VrV HAr 
o z, o o o 2*i e z ’ e e e 


iH = H o- H e = — 


Z r V Ar y r V_ . Ar 
o z, o o e z,e e 


(B6) 


Mass-averaged energy addition (ideal head rise): 

;E r o V z, o U o V 0 , o i t I e V *, e U e V e • « ir e 


a _ U o V 6,o U e V 0,e . gr h 
1 g g r t 


h 


E r V Ar 
o z, o o 


V r„V„ „ Ar, 


e z, e e 


h 


(B7) 
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Mass-averaged efficiency: 



Mass-averaged head coefficient: 




gAH 



Average flow coefficient (no boundary layer correction): 


(B8) 


(B9) 


1 = 


Qy 

448. 8 tr(r 2 - r 2 )U t> x 


Cavitation number: 


A local cavitation number is defined by 


_ h e - ly _ 2 gH»v v e 

(v ^) 2 (v ;) 2 (vp 2 

2g 


(BIO) 


(B12a) 


For application to the inducer where the inlet total head is constant radially and the inlet 
flow has no whirl component (Vg = 0) an average or blade cavitation number of the 
following form is used: 

f: _ 2 S H sv »e 

u ?,e< 1+ ^> 1+ ^e 


(B12b) 
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